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ION, ELECTRON AND NEUTRAL TEMPERATURES DERIVED FROM 
ION COMPOSITION DISTRIBUTION 


B. C. Narasinga Rao* 
Laboratory for Space Sciences 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 


Recently Brinton et al (1968) reported the positive ion 
mass spectrometer measurements obtained from the Geoprobe 
(Argo D-4 rocket) over an altitude range of 200-630 km at 
1300 EST on March 2, 1966 above Wallops Island, Virginia. 

The ionic constituents measured were 0 + , N + , H + , He + N0 + , 

0g and Ng (Fig. 1). The measurements showed 0 + as the major 
constituent throughout the altitude range. Since the dis- 
tribution of each ionic species is controlled by diffusive 
processes at higher heights and chemical processes at lower 
heights with an overlap in the intermediate heights (Bauer 
1964) it is possible to derive the ion, electron and neutral 
temperatures from the composition data using the appropriate 
expressions . 

Under diffusive equilibrium, the distribution of an 
ionic species X + in an ion mixture having a mean ionic mass 
m + is given by (Mange 1960, Bauer 1966) 
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where T ~ electron temperature 
0 

Tj. ion temperature (assuming all ions to be at the same 
temperature) 

z 1 - geopotential height 

g Q •« acceleration of gravity at the earth’s surface 
As pointed out by Bauer (1964, 1966) an important feature 
of the diffusive equilibrium for ions is that, whereas the 
distribution of an individual ion species depends on T as 
well as 1\ , the ratio of the two species in diffusive equilibrium 
depends on T\ only and, thus, the relative abundance of ionic 
species in diffusive equilibrium provides a measure of T^. 

From eq. (1) we obtain 


n(X*) m(Xo)-m(xJ 

v In [ i-] - 2 i- (2) 

n(X+) WT ± /g 


The above equation can be used to derive from the height 
gradient of the ratio of n(X^) and n(Xg), if both of them 
are in diffusive equilibrium. If they have not attained complete 
diffusive equilibrium in the region under study, the temperature 
derived from such a profile will be in error depending on the 
extent to which the distribution deviates from that of 
diffusive equilibrium. From the nature of distributions in 
Fig. 1 it appears that 0 + , H + , N + and He + are probably in 
diffusive equilibrium from 500 km upwards. From theoretical 
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studies it is estimated that whereas 0 + , H + and He + attain 
diffusive equilibrium by about 500 km or even lower altitude 
(Nisbet 1963, Bauer 1964, 1966; McElroy, 1965), N + is still 
under chemical-cum-dif f usive control (Bauer, 1966). It is 

*4“ -j- 

thus reasonable to proceed with the assumption that 0 , H 
and He + are in diffusive equilibrium above 500 km to determine 
using eq. (2). In fact, even for n(H + ) the theoretical 
scale heights calculated on the basis of diffusive equilibrium 
and using the temperatures and their gradients presented 
later, are 130 km, 180 km and 320 km which are close to the 
observed values of 100 km, 170 km and 320 km at altitudes of 
500 km, 550 km and 600 km respectively. 

The values of derived from the height distributions 
of n(H + )/n(0 + ) and n(He + )/n(0 + ) using eq 0 (2) are shown in Fig. 
2 for the height range of 350 to 600 km. It is seen that from 
about 450 km the T i values derived from the two ratios are 
essentially the same. The mean values are 1800°K at 600 km, 
1500°K at 550 km, 1200°K at 500 km and 900°K at 450 km. 

These values are comparable to the T A values obtained by 
Thompson scatter experiment by Evans (1967) under similar 
conditions. The T^ values below 450 km derived from 
n(He + )/n(0 + ) distribution are in the range of 250°-350°K 
and are much too low compared to neutral temperature, though 
the T^ values derived from n(H + )/n(0 + ) are comparable to 
neutral temperature. This result suggests that the n(He + ) 
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distribution is not in diffusive equilibrium below 450 km 
while n(0 + ) distribution may correspond to diffusive equilibrium. 
The n(H H ) distribution, even if it is under chemical equilibrium, 
may still look like that of diffusive equilibrium, since both 
the distributions are similar for n(H + ) when it is a minor 
ion (Bauer 1966a) . 

We may utilize these values of and its height 
variation to derive T 0 as follows. The geopotential scale 
height H* of the electron density distribution under diffusive 
equilibrium is given by 


H* 


T +T . 
e i 
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Since H* and m + can be determined from the composition data, 

eq. (3) may be used to derive the T distribution. The calculation 

of T is performed by a method of successive iteration starting 

from ?)T /?s z* = 0. The T values thus obtained are shown in 

e e 

Fig. 2. They are 2900°K at 600 km, 2870°K at 550 km, 2790°K 
at 500 km and 2580°K at 450 km. These values are in close 
agreement with the measured values of T by the electrostatic 
probe experiment by Brace et al (19G8) on the same flight. 

These experimentally measured values are also shown in Fig. 2 
for comparison. 
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With the values of T 1# T^, STj/az* and rfl^/dz*, thus 
obtainedi it is possible to determine if the n(N + ) distribution 
is in diffusive equilibrium at higher altitudes. It is 
found that at 600 km, the scale height of the observed dis- 
tribution is 275 km whereas the diffusive equilibrium scale 
height is 185 km. This difference implies the lack of 
diffusive equilibrium. However the distribution of n(N + ) 
can be used to calculate T n . Bauer (1966) discussed the 
problem of n(N + ) distribution below and above its peak and 
showed that, above the peak, it is controlled by production 
through charge transfer between He + and Ng and loss through 
the divergence of diffusive flux, and at still higher heights, 
by diffusive equilibrium. Above the peak and below the 
equilibrium region it is approximately given by 

n(N + ) « i n(He) « <sxp(- ■§') (4) 

where K (a parameter dependent on diffusion velocity) may 
be assumed constant. Then the n(N + ) distribution follows the 
neutral helium distribution, which in turn, depends only on the 
neutral temperature. The values of T n derived from the n(N + ) 
distribution in the height range of 350-450 km is 760°K. 

It is also possible to determine T from the nCHe" 4 ”) 

n 

distribution in the height range of 350-400 km where it is 
mainly controlled by chemical processes. Under chemical 
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equilibrium n(He h ) is given by (Bauer 1966) 

. I n (He) , 

n(He ) — oc exp(wj) (5) 

Y n(N 2 ) ^ 

where I is the photoionization rate coefficient and v is 
the charge exchange rate coefficient of He + with N 2 . In the 
above equation the charge exchange reaction of He with Og for 
loss of He + is not considered since n(0 2 ) is lower than n(Ng) 
by an order of magnitude in this altitude range. Then n(He + ) 
increases with height with a scale height corresponding to the 
neutral temperature and mass number 24. The value of T n de- 
rived from n(He + ) distribution in the height range of 350-400 km 
is 780°K. This value is very close to that derived earlier 
from the n(N + ) distribution and is also close to the value of 

T derived from neutral N« distribution as measured on the same 
n 2 

flight by Pelz and Newton (1968) . 

In conclusion it is noted that the ion, electron and 
neutral temperatures may be derived from the ion composition 
data in a consistent way. 
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CAPTIONS FOE FIGURES 


Fig. 1 


Fig . 2 


Ion density distributions observed over Wallops Island , 
Virginia by positive ion mass spectrometer aboard the 
Geoprobe (Brinton et al, 1968). 

Ion, electron and neutral temperatures derived from ion 
composition data show., in Fig. 1. The T 0 and T R 
distributions shown by continuous lines are obtained 
by electrostatic probe by Brace et al (1968) and derived 
from neutral N 2 distribution measurements by Pelz and 
Newton (1968) respectively made on the same flight. 


WALLOPS ISLAND, VA. MARCH 2, 1966 1800 UT 1300 EST 
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Figure 2 
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